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Abstract. New additional equations for the Newtonian dynamical systems on Riemannian manifolds are found. They 
supplement the previously found weak normality conditions up to the complete normality conditions for Newtonian dynamical 

■ systems. 

Si 

^ ■ 1. Introduction. 

Let M be the Riemannian manifold with metric tensor <?jj . The Newtonian dynamical system on M is the law of 

■ motion of particles on M given by the following equations in local coordinates 

o\ : 

■ (1.1) ± l =v l = &(x 1 ,...,x n ,v 1 ,...,v n ) 

o 

The equations (1.1) realize the Newton's second law for the mass point with unit mass. The concept of normality for 
Newtonian dynamical systems was introduced and then investigated in the series of papers [1 - 5] . It is based on the 
, study of some hypersurface S in M. For some point P on S we define the initial data 

^ (1.2) x% =Q =x\P) v% =o = v{P)n\P) 

for the equations (1.1). Here n l {P) are the components of unit normal vector at the point P and v(P) in (1.2) is some 
scalar function on S. The equations (1.1) and initial data (1.2) determine the particle flow starting at t = from S 
along the normal vector n(P) with the initial velocity v(P). For t > these particles form another hypersurface St 
L* ' and determine the family of diffcomorphisms f t : S — > S t . 



Definition 1. The family of diffeomorphisms f t is called the normal shift along the trajectories of dynamical system 
if the trajectories of (1.1) cross all hypersurfaces St along their normal vectors. 

Definition 2. Newtonian dynamical system (1.1) is called the dynamical system accepting the normal shift of hyper- 
surfaces if for any hypersurface S C M one can find the function v(P) on S such that the family of diffeomorphisms 
ft given by (1.1) and (1.2) is the normal shift of S. 

The normality condition from definition 2 was first stated in [1] and [2]. By the analysis of this condition in [3] for 
Euclidean case M = M" two relatively independent normality conditions were derived: weak normality condition and 
additional normality condition. Each of them is written in the form of the system of partial differential equations for 
the functions <& l from (1.1). They two both form the sufficient condition for the normality condition from definition 
2 to be fulfilled for the dynamical system (1.1). 

In [5] the weak normality condition was generalized for the non-Euclidean case of an arbitrary Riemannian manifold 
M. In this paper we generalize the additional normality condition from [3] for the case of Riemannian manifold and 
we get the complete normality condition for (1.1) in form of aggregate of weak and additional normality conditions. 
Here as in the Euclidean case of W 1 the complete normality is sufficient but not necessary for the normality condition 
from the definition 2. 

Since it seems difficult to obtain the conditions equivalent to definition 2 and written in form of differential equations 
for the functions <I> 1 we shall replace the definition 2 by more strict definition 3. 
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Definition 3. We say that dynamical system (1.1) satisfies the strong normality condition if for any hypersurface 
S C M, for any choice of the point P e S and for any real number w 7^ one can find the function v(P) such 
that it is normalized by v(P ) — v and the family of diffeomorphisms defined by this function is the normal shift of 
hypersurface S. 

Note here that for to avoid the multivalued functions v(P) one should always consider only connected and simply 
connected hypersurfaces in definitions 2 and 3. 

2. ON THE EXPANSION OF THE ALGEBRA OF TENSOR FIELDS. 

Systems of differential equations of the form (1.1) are usually connected with the vector fields on the manifolds. In 
our case the right hand sides of (1.1) form the vector field on the tangent bundle TM for the manifold M. Digressing 
for a while from the particular vector field given by (1.1) let us consider some vector V tangent to the tangent bundle 

V = X 1 — H hi — + W 1 —— H h W n — 

dx 1 dx n dv 1 dx n 

First n components of the vector (2.1) form the tangent vector X = 7r(V) to M. Other components of (2.1) after 
some modification form another tangent vector Z = p(V) to M whose components are Z l = W % + Y\ jk v k XK Two 
linear maps 7r and p project the vector V onto the pair of vectors X and Z tangent to the manifold M. 

Projections 7r and p applied to the vector fields on TM however will not give the vector fields on M. They give 
the elements of quite other set: the expanded algebra of tensor fields on M. Exact definition of tensor field from such 
algebra is the following. 

Definition 4. Tensor field of expanded algebra is a map argument of which is a point of tangent bundle TM and 
the value of which is in the tensor algebra build over the tangent space to M at the point being the projection of 
argument from TM to M. 

Expanded algebra of tensor fields is equipped with the natural operations of tensor product and contraction. It is 
also equipped with two operations of covariant differentiation Vj and Vj . The detailed discussions of these features 
of expanded algebra of tensor fields can be found in [5]. 

Each point of TM is the pair of point P G M and tangent vector v at this point (vector of velocity). If we map 
each point of TM onto the corresponding vector v then we get the vector field v of expanded algebra. From v we 
construct the scalar field v = |v| of expanded algebra being the modulus of velocity. In addition to these two fields we 
define the vector field N = |v| _1 v that consists of unit vectors directed along the vector of velocity v. It is defined 
only at that points of TM where v ^ 0. By means of the components of N we construct two projector valued fields 

(2.2) Ql = N k N* PI = (SI - N k N*) 

from the expanded algebra. Various relationships with projectors Q and P from (2.2) can be found in [5]. 

Replacing time derivatives of velocity in (1.1) by its covariant time derivatives we can rewrite (1.1) as follows 

(2.3) x i = v i V t v i = F i (x 1 ,...,x n ,v 1 ,...,v n ) 

where F l = +Tj k v k v : > are the components of vector field F of expanded algebra known as a force field of Newtonian 
dynamical system (2.3). 

3. Weak and complete normality conditions. 

For the dynamical system (2.3) we consider the Cauchy problem with initial data (1.2) on some hypersurface S. 
Let us choose the local coordinates u , . . . 1 u n ~ 1 on S. Solving the above Cauchy problem we obtain the family of 
trajectories of the dynamical system (2.3) on M parameterized by u 1 , . . . , u 11 ^ 1 . Variation of these variables defines 
the following vectors n tangent to M 

dx 1 d dx 11 d 

du* dx 1 du 1 dx n 

on the trajectories of (2.3). Scalar products of T\, . . . , r„_i and N are the scalar functions ipi = (rj,N) = gk q r k N q . 
These functions define the mutual orientation of trajectories and the hypersurfaces St- 



(2.1) 



(3.1) 



COMPLETE NORMALITY CONDITIONS 



3 



Definition 5. Say that Newtonian dynamical system (2.2) satisfies the weak normality condition if each function 
ipi = tpi(t) is the solution of linear homogeneous ordinary differential equation of second order for any choice of 
parametric family of trajectories of it. 

The main result of [5] is that the weak normality condition is equivalent to the following system of nonlinear partial 
differential equations for the force field F of dynamical system 

(v- 1 F l + WdF k N k ))Pl = 

(3.2) (V,F fe + V fe F, - 2v- 2 FiF k )N k pi+ 

+ v-^VuF^ - V k F r N k N r FJI» = 

Let the weak normality condition in form of the equations (3.2) be fulfilled. Then for to get normality in the sense of 
definitions 2 and 3 we should provide the following initial data 

(3.3) ¥>4=o = ° H=o = ° 

for the functions ipi by the proper choice of function v(P) = . . . ,u n_1 ) in (1.2). There n(P) is the unit normal 
vector for the hypersurface S. From (1.2) we obtain 

(3.4) Nl^^u 1 ,...,!*"- 1 ) 

Because of (3.4) first part of the initial conditions (3.3) is satisfied for any choice of v(P). Now let us proceed with 
second part of initial conditions in (3.3) 

(3.5) &| t=0 - V t rf ; Nj\ t=0 +g jk T?V t W\ t=0 

For VtNi we use (3.15) and (3.19) from [5] and then we obtain V t 7V J = v~ 1 P^F q . For the covariant derivatives of 
the vectors (3.1) we use (3.18) from [5] 

<9V . _, dxP a . 



lt=0 dtdu llt=0 pq du l 14=0 
Taking into account the initial data (1.2) we may transform this relationship into the following form 

(3-6) v 4 t/U = ^| t=0 + ^| t=0 + r>V | t=0 

For the further transformations of the equations (3.6) we should recall some facts concerning submanifolds of Ric- 
mannian spaces 

(3-7) ^ + W = r? k 1 + W 

(3-8) %+T> q -! = -*H 

Here bik and bf are components of tensor of second quadratic form and T q ik are components of metric connection on S 
defined by the metric <jik = gpqrf r| induced from M to the hypersurface S. The equations (3.7) and (3.8) are known 
as Gauss and Weingarten formulae (see [6] and [7]). Comparing (3.6) with (3.8) we get 

(3-9) V t Tf| t=0 = ^| t=0 -^| t=0 

Now we substitute the above obtained formula Vt-ZW = v~ 1 P^F q and (3.9) into (3.5). Then from (3.3) we obtain the 
following equation for still unknown function v — ^(w 1 , . . . , w" -1 ) 

(3-10) ^ = -tr 



4 BOLDIN A.YU., BRONNIKOV A.A., DMITRIEVA V.V., SHARIPOV R.A. 

Components of the force field in (3.10) depend on the velocity vector therefore they contain the dependence on the 
unknown function v in the form of v p = ^(m 1 , . . . , u n ~ 1 )n p (u 1 , . . . , u n_1 ). 

The equations (3.10) form the overdetermined system of differential equations. For to have common solution v they 
should satisfy the compatibility conditions. We obtain these conditions when consider the following derivatives 

calculated according to the equations (3.10). To make shorter all further calculations we define covariant derivatives 
Di as covariant derivatives with respect to u l given by the formula (3.18) from [5]. These covariant derivatives are 
applicable to the tensor-valued functions whose domain of definition is S. For the tensor fields of M restricted to the 
hypersurface S these covariant derivatives are calculated as Di = T k V k . Note that Di are not applicable to tensor 
fields of expanded algebra unless some lifting of S from M to tangent bundle TM is defined. For the force field F in 
(3.10) and (3.11) such lifting is given by the second part of initial conditions in (1.2). Therefore 

(3.12) DiF p = r k V k F p + {n k D lV + vD in k )V k F p 

The results of applying Di to vector fields T 3 k and n J are defined by Gauss and Weingarten formulae (3.7) and (3.8) 

(3.13) Di4 = f ^ + b ik ni Div? = -b\r{ 

So the derivatives Di establish the link between inner geometry of S and the geometry of outer space M itself. 

The derivatives DiV for the scalar function v on S are given by the equations (3.10) the compatibility condition for 
which we are going to derive now. Let's combine (3.12) and (3.13) 

' = -v- 3 F v TfF a T? - v- l V v F a rfrl- 



(3.14) dv>dv?~ P i Q j ^I'i'i 

+ v- 2 F pT f V r F q n r r] + V p ^6[r r p r? - v^T^Fg - v-%n^F q 

After exchanging indices i and j in (3.14) we obtain another expression for the same derivative in left hand side of 
(3.14). The difference of these two expressions should be zero. This gives us the compatibility condition for (3.10) in 
the following form 



(3.15) 



rM f V g F fc - VfcF, + Nr V r FgF k -V r F k F q ^ + 



v 

+ v:T k r Tp k F q -V jT i;T?V k F q = {) 
In order to simplify the equations (3.15) we need to recall the following relationships due to (3.4) 
(3-16) {g krT l^) T 1=Pl 

By means of contracting (3.15) with the quantities enclosed in brackets in (3.16) we obtain 

pk p q ( Vg-Pfc ~ ^kF q + Nr V r F q F k - V r F k F q \ + 



(3.17) 



v 



+ H k P]V k F q - H k P?V k F q = 



where H k are the components of the symmetric operator H determined by the second quadratic form of hypersurface 
S 

(3-18) H* = g lr r r q ^bP J T k p 

Symmetric operator H with the matrix (3.18) satisfies the following relationships 



(3.19) 



HP = PH = H 



rank(H) < n - 1 
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Lemma 1. For any point P in M and any operator H in tangent space to M at this point satisfying the properties 

(3.19) one can find the hypersurface S passing trough the point P such that the matrix of the operator H is given by 
formula (3.18) at this point. 

Proof of this lemma technical. We omit it noting only that one can choose S in the class of quadrics for some local 
coordinates x 1 , . . . , x n in the neighborhood of P. 

Lemma 1 shows that the matrix of the operator H in (3.17) is rather unrestricted. This let us replace (3.17) by 
two separate equations of the following form 

P?PJ ( V q F k V k F q + N r VrF q F k -V r F k F q \ = Q 

(3.20) J \ v J 

H*P]V k F q - H*P?V k F q = 

Now consider the operator K with the matrix K\ = P^V ' k F q P^g r K Its properties are similar to (3.19) i.e. 

(3.21) KP = PK = K rank(K) < n - 1 

Second equation (3.20) then means that the product KH is symmetric operator. Because of lemma 1 the operator H 
is arbitrary symmetric operator. Taking H = P and using (3.21) we get that K is also the symmetric operator. The 
product of two symmetric operators KH is symmetric if and only if they are commutating. Thus we have KH = HK 
for arbitrary symmetric operator H satisfying the relationships (3.19). This is possible only if K is proportional to P 
with some scalar factor K = AP. Scalar factor A is easily calculated as A = tr(K)/tr(P). Now we can write 

(3.22) V / V / 

P k X7, F q P r 

P ^ kFqplq = ^i k Y" pf 

excluding the matrix H from (3.20) at all. The equations (3.22) just above combined with (3.2) form the complete 
normality conditions which are sufficient for the definitions 2 and 3 to be fulfilled for the dynamical system (2.2). 
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